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ABSTRACT 

Aims. To explore the solar coronal heating enigma by an analytical and numerical study of the enhanced phase mixing of harmonic 
Alfven waves propagating in gravitationally stratified coronal structures of varying magnetic field divergence. 

Methods. Corrected analytical solutions are derived to model the dissipation of Alfven waves propagating in divergent and stratified 
coronal structures. These analytical solutions are validated and further explored using a newly developed 2.5D visco-resistive linear 
MHD code. 

Results. Corrected analytical solutions describing the enhanced phase mixing of Alfven waves in divergent and stratified coronal 
structures are presented. These show that the enhanced phase mixing mechanism can dissipate Alfven waves at heights less than half 
that is predicted by the previous analytical solutions. In divergent and stratified coronal structures, the enhanced phase mixing effect 
occurs only when the ratio of the magnetic and density scale heights, H b /H p < 2. The enhanced phase mixing of 0.1 Hz harmonic 
Alfven waves propagating in strongly divergent, H\, = 5 Mm, stratified coronal structures, H p = 50 Mm, can fulfill 100% of an active 
region heating requirement, by generating viscous heating fluxes of Fu as 2100 J nr 2 s -1 . The Alfven waves in this configuration are 
fully dissipated within 20 Mm, which is six times lower than would occur as a result of standard phase mixing in uniform magnetic 
fields. This results in the heating length scale, Lh, defined as the height at which 95% of the Alfven wave poynting flux has dissipated, 
being lowered by a factor of six, to less than half of an active region density scale height. Using the corrected analytical solutions it 
was found that, for a given wave frequency, the generation of a heating length scale of Lh < 50 Mm, by enhanced phase mixing in 
strongly divergent magnetic fields, requires a shear viscosity eight orders of magnitude lower, than required by standard phase mixing 
in uniform magnetic fields. It was also found that the enhanced phase mixing of observable, co as 0.01 rads s , Alfven waves, in 
strongly divergent magnetic fields, H b = 5 Mm, can generate heating length scales within a density scale height, H p = 50 Mm, using 
classical Braginskii viscosity. It is therefore not necessary to invoke anomalous viscosity in corona, if phase mixing takes place in 
strongly divergent magnetic fields. This study shows that the importance of enhanced phase mixing as a mechanism for dissipating 
Alfven waves in the solar corona (a stratified and divergent medium), has been seriously underestimated. 
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1. Introduction 

Phase mixing was proposed by IHevvaerts & Pries! d 1983b as a 
mechanism for dissipating Alfven waves in the solar corona. 
They suggested that shear Alfven waves propagating on neigh- 
boring magnetic field lines, in inhomogeneous plasmas, would 
move out of phase with each other. This would then lead to 
strong gradients, or short length scales, developing perpendic- 
ular to the direction of wave propagation. The build up of strong 
transverse gradients eventually leads to dissipation of the wave's 
energy via shear viscosity or resistivity, thus heating the plasma. 

The ability of phase mixing to make a significant contribu- 
tion to coronal heating has attracted ample attention of solar 
physicists. This simple mechanism has been applied to a variety 
of coronal structures; most notably coronal holes, loops and ar- 
cades (see reviews in lBrownind[l99U iNara in & Ulmschneider 
119901 119961: iParkerl |1991|). These studies are supported by 
observational evidence of Alfven waves, which after being 
generated by photospheric footprint movements, p ropagate up 
through these coronal structures (see llrelandll 1 9961) . Phase mix- 
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ing may also contribute to coronal heating through the pro- 
duction of fast and slow magnetosonic waves due to non- 

i ii 1 

linear coupling with Alf ven wav es (see Nakariakov et al. 1997; 

iTsiklauri & Nakariakovl 120021: iTsiklauri et al.1 120021 120031 
Recentlv lTsiklaurild2006al) has shown that this mechanism could 
lead to significant coronal heating through the generation of par- 
allel electric fields. The main alternative to coronal heating by 
waves, co mes from the well studied reconnec tion of magnetic 
fields (see lPriesdl2TjQ3r.|Priest et al.ll2003ll2005h . 

While Alfven waves contain enough energy to account for 
the observed one million kelvin temperatures of the corona, 
their notoriously slow dissipation means that their energy is car- 
ried high into the corona. The low classical Braginskiil (1 1965b 
values of shear viscosity and resistivity mean that, even with 
phase mixing, it is very difficult to dissipate Alfven waves 
within the few density scale heights, H p « 50 - 200 Mm, re- 
quired for them to be responsible for coronal heating. Numerous 
studies have therefore sought to lower the height at which 
Alfven waves dissipate by a variety of mechanisms: localizing 
the Alfven w aves into closed-field reg ions (e.g. Resonant ab- 
sorption (see lOfman et al.1 11994 1 19951) ). using large (anoma- 
lous) values of resistivity and viscosity, and by what we will 
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call the enhanced phase mixing effect. In this latter mecha- 
nism, the phase mixing of harmonic Alfven waves propagat- 
ing in divergent magnetic fields was foun d to generate dissi- 
pation rates far in excess of the standard iHeyvaerts & Priest! 
(1983)) exp(-z 3 ) rate (where z is height). Numerous authors 
have therefore investigated how magnetic field divergence, along 
with gravitational density stratification, cha nges the efficiency 
of phase mixing to di s sipate Alfven waves (Poedts et al.lll997t 
IRuderman et al.|ll998l [19991: iDe Moortel et al.l 1 1 999l l2000h . A 
recent study by lAschwanden et al.l (120071) suggests that 83% of 
the c orona's heating requireme nt is located in active regions, 
while Aschw anden et al. d2000l) showed that the heating scale 
height is required to be less than the active region's density scale 
height. Therefore, for Alfven waves to be responsible for the ma- 
jority of coronal heating, a mechanism that can dissipate them in 
act ive regions, within a den sity scale height, is required. 

IDe Moortel et alj d2000t) detail the effects of altering the den- 
sity scale height on the phase mixing of Alfven waves prop- 
agating in radially divergent magnetic fields. They concluded 
that the resultant dissipation can either be enhanced or di- 
minished, depending on the specific configuration of the coro- 
nal structure. Analytical solutions for Alfven waves propagat- 
ing in more gene ral divergent and stratifi ed coronal structures 
can be found in IRuderman et al.l dl998l) . In non-stratified di- 
vergent coronal structures, they found that the standard dis- 
sipation rate of harmonic waves, exp(-z 3 ), is replaced by an 
even faster exp(- exp(z///b)) rate which depends on the mag- 
netic scale height, H\,. In stratified coronal structures, harmonic 
Alfven waves were generally found to dissipate slower than the 
standard rate; although this depended strongly on the specific 
coronal structure. These two papers used quite different models, 
making a direct comparison difficult, but they do both generally 
agree that a diverging magnetic field enhances the phase mixing 
mechanism, while the gravitational density stratification dimin- 
ishes it. They also raise a number of important questions: what 
are the coronal conditions necessary for enhanced phase mixing 
to occur? How much lower do harmonic Alfven waves dissipate 
in the corona as a r esult of enhanced phase mixing? Can the ana- 
lytical solutions of lRuderman et al.l (Il998h be confirmed numer- 
ically and what is their range of validity? Finally, is enhanced 
phase mixing a viable mechanism for heating the solar corona to 
temperatures in excess of one million kelvin? 

The aim of this study is to attempt to answer these questions. 
In the next section we present the model used in our study. In 
Sect. 3 we correct the analytical solutions of Ruderman et alJ 
(1998), which describe the phase mixing of Alfven waves in 
stratified coronal structures with divergent magnetic fields in the 
WKB approximation. In Sect. 4 we detail the numerical code 
that was used to solve the linearized MHD equations which de- 
scribe phase mixing. In Sect. 5 we present our numerical results 
and compare them to the analytical solutions of Sect. 3. Finally 
in Sect. 6 we discuss the results and important conclusions that 
can be drawn from this study. 



2. Analytical Model 

The starting point for our analysis is the system of MHD equa- 
tions for a cold, incompressible plasma 

dV 1 

p— + p(V ■ V) V = B x(VxB)-Vx (pvV x V) , (1) 

at fi 



Here B is the magnetic field, V the plasma velocity, p the mass 
density, v the kinematic viscosity, rj the magnetic diffusivity and 

the magnetic permeability of free space. In the momentum 
equation, only the shear viscosity term is included, even though 
the bulk viscosity in the solar corona is many orders of mag- 
nitude larger. This approximation is valid as the bulk viscos- 
ity is related to compressibility, which is not included in our 
model and thus does not effect the amplitude of propagating 
Alfven wa ves. Indeed, it has be en s hown through numerical sim - 
ulations bv lOfman etail d 19941) and lErdelvi & Goossensl d!995l) . 
that, under typical coronal conditions, the bulk viscosity can be 
neglected in Alfven dissipative layers in comparison with the 
shear viscosity. The use of the cold plasma approximation in our 
model applies only to the perturbation quantities of velocity and 
magnetic field. This approximation holds in the solar corona, as 
the plasma beta is usually very small (/3 « 0.01). The equilibrium 
quantities incorporate thermal physics to account for the vertical 
density stratification of the corona. 

The applicability of the MHD equations for the description 
of wave propagation in the solar corona, is very often questioned 
on the grounds that the coronal plasma is collisionless. In this 
respect, we have to note that it is meaningless to discuss if a 
plasma is collisional or collisionless before we have specified 
the characteristic time scale of the problem. A plasma has to 
be considered as collisionless if the characteristic time scale is 
much smaller than the ion collision time. Conversely it has to be 
considered as collisional if the characteristic time scale is much 
larger than the ion collision time. The collisio nal time of pro - 
tons in the solar corona, r p , was calculated bv lHollweil 0985). 
He obtained that t p can vary from about 0.7 s in active regions 
to about 7.5 s near the base of coronal holes. For the partic- 
ular values of the electron number density (6.43 x 10 14 irT 3 ) 
and temperature (2 x 10 6 K) used in our numerical calculations 
we obtain r p s 3 s. In the following calculations, we consider 
waves with periods larger than or equal to 10 s, so that the use 
of the MHD description is fully justified. It is also known that 
the relatively strong magnetic fields present in the corona, in- 
hibit cross-field particle motions, while wave particle interac- 
tions are known to impede particle motions even along the field. 
These have the combined effect of localizing the particle inter- 
actions in all directions, thereby enabling a fluid description to 
apply (see iPriest & Forbes! 120001) . Given this, there have been 
various attempts to include e ffects beyond MHD into coronal 
studies; Tsiklauri et al. (2005) included ion and electron kinetic 
effects, Ofman & Vinas (2007) considered multiple species of 
ions, while iTsiklauril d2007l) has included the effects of two fluid 
species. 

In what follows we use the Cartesian coordinates x,y,z to 
represent respectively the transverse coordinate, the ignorable 
coordinate {d/dy = 0) and the height in the corona. We neglect 
the magnetic field diffusion (77 = 0). We assume a static (Vq - 0) 
two dimensional (Z?oy = 0) equilibrium, where the subscript '0' 
indicates an equilibrium quantity. We then linearize Eqs. (fl}- 
d2) and consider perturbations in the form of linearly polarized 
Alfven waves, so that only the y-components of the perturba- 
tions of V and B are non-zero. This results in the following sys- 
tem of linear equations: 



Po- 



dVy 

~di~ 



1 

Mo 



(Bo-V)B y + V-(povW y ), 



(3) 



— = Vx(VxB)-Vx (77V x B) . 
ot 



(2) — £ = CB -V)V y . 
ot 



(4) 
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Eliminating B y from Eqs. (O-© gives the diffusive wave equa- 
tion 



d 2 v v 



Po- 



-- (B ■ V) 2 V y + V ■ LvV^ 



(5) 



Since the equilibrium magnetic field is two-dimensional and 
divergence-free, it can be expressed in terms of the flux function 



dip dip 
Box = —Boo-z-, B 0z = Boo—, 
dz dx 



(6) 



where Z?oo is the magnitude of the magnetic field at the coor- 
dinate origin (x = z = 0). An individual magnetic field line is 
described by ip — const. We can introduce a second function, (p, 
satisfying Vip • V(p = 0, and consider ip and <p as new curvilinear 
coordinates in the xz-plane. Thus Eq. (0 can be re-written as 



(T- 



8 2 Vy 



v z A0 J 
+4 

dt 



ddVy 

d(f> d(f> 



dip\ * dip j d<p\ * dip 



(7) 



where cr - po/poo is the dimensionless density, poo is the density 
at the coordinate origin, Vao is the Alfven speed at the coordinate 
origin given by Vao = Boo /(/Woo) > J is the Jacobian of the 
coordinate transformation, 



dip dtp dip dip 
dx dz dz dx' 

and hj, and h$ are the scale factors given by 



\diff) + W 



1/2 



dz \ I dx 

a* + 



1/2 



(8) 



(9) 



N ote that our defin i tion o f these scale factors differs from those 
of Ruder man et al.l d 1998b . resultin g in a slightly different form 
of Eq. . The scale factors used by Ruderman et al. ( 1998), 
and h$r, are related to ours by = Jh^ and - Jh^. In either 
case, both definitions of the scale factors lead to identical general 
solutions. 



3. Analytical Solutions 

iRuderman et al.l d 19981) obtained an analytical solution of Eq. 
in the WKB approximation. Their approach is valid when the 
following assumptions are satisfied 

1 . The ratio of the characteristic scales in x and z-directions is 
small, xo/H <k 1; 

2. The wavelength is < xq\ 

3. The characteristic scale of damping is < H. 

Here H is the smaller of the two quantities, the characteristic 
scale of the magnetic field variation, H^,, and the density scale 
height H p . 

We verified the calculations of IRuderman et al] d 1998b . and 
found that they are valid up to their Eq. (36). However there is 
an error in their Eq. (37). The correct equation should be 



I(x,z) 



z dx 1 dVp 



dipvl 



dx 



dz', 



(10) 



where the Alfven speed is given by 
Bo 

V A 



(ID 



0"oPo) 1/2 ' 

Ruderman et al.l (1 19981) missed the factor dx/dip in their expres- 
sion for I(x, z). If we take V y oc exp(-/wf), then the evolution of 
the wave velocity is given by 



V y (x,z) = o-- l/4 W(x)e 



,-A(«) 



where 



Mx,z) 



2B 00 



r 

Jo 



v(x,z')B 2 Jx,z')I 2 (x,z') 



V A (x,z') 



dz', 



(12) 



(13) 



and W(x) is determined by the boundary condition at z = 0. For 
a detailed derivation of this general solution, see Appendix A. 

Now, using the corrected expression for I(x,z') we re- 
consider three particu lar examples of equilibria studied by 
Rude rman et al.l JT998). Following them we assume that v = 
const. The first equilibrium they studied was a uniform verti- 
cal magnetic field in an isothermal atmosphere, so that Z?o x = 0, 
Z?0z = Boo and po(x,z) = po(x)e~ z/Hp . Where again, the term 
isothermal applies to the equilibrium quantities only. For that 
specific configuration iff = x, <p = z, dx/dip — 1 and the refore 
the expression for I(x, z') given by IRuderman et all ((1998) coin- 
cides with Eq. ( [TOl l. Hence, in this particular case, the analytical 
solution obtained by Ruderman et al. is correct. 

The two other cases studied by IRuderman et all (Q998) were 
equilibria with (i) constant density along z (non-stratified) and an 
exponentially diverging magnetic field, and (ii) constant Alfven 
speed along z- We consider a more general equilibrium consist- 
ing of an exponentially diverging magnetic field in an isothermal 
stratified atmosphere. In this case the magnetic flux function is 
given by 



ip(x,z) = H b e- Z/Hb sm(x/H b ), 



(14) 



so that the magnetic field is determined by 

Box = B m e- Z/Hb sin (x/H b ) , B 0z = B aa P zlH <' cos (x/H h ) . (15) 



(16) 



The orthogonal curvilinear coordinate <p is given by 
(p(x,z) = H h e-' JHb cos(x/// b ). 
The expression for the equilibrium density is 

p (x,z)=Mx)e- z/H *. (17) 

The exact form of the missing factor from Eq. ( fTOb is given by 
dxfdifi = sec (x/Hb) exp(z/H\,). In the following solutions, we 
use the simplified form dx/dip ~ exp(z/Hb), since sec (x/H\,) * 
1 . Then, using Eq. ( fTOl ), we obtain 



I(x,z) = 



H p H b [l - exp(2z/// b - z/2H p )] d p Q 
^0P /2 p ^ 2 (4i/ p -//b) dx ' 
Substituting Eq. ( [TBI into Eq. ( IT3l we arrive at 



(18) 



A(x,z) » A(x) 



1 - exp(3z/# b - 3z/2# p ) 



3 (ff b - 2H P ) 
l-exp(z/H h -z/H p ) 

Hp - Z/b 
1 - exp(-z/2ff p - z/gb) 
Hb + 2H p 



(19) 
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A(x) = 



VLO 



T/ 3 .1/2 3/2 , 

y AO Po Poo (4H P 



\2 



dpo 
dx 



(20) 



Clearly, Eqs. C[9]l-(|20]) are only valid when H h + H p ,H b + 2H p 
and H b + 4H p . In what follows we assume that these inequalities 
hold. 

By considering a non-stratified coronal structure permeated 
by a uniform magnetic field, where H a — > oo and H b — > oo, we 
can obtain the Hevvaerts & Priesil d 19831) standard phase mixing 
solution; Eqs. (fT9b — d20l> reduce to 



A(x,z) 



6V 



dV A 
dx 



(21) 



To obtain the equilibrium which represents a non-stratified 
divergent coronal structure, we take H p — > oo. Then Eqs. (fT9l- 
( l20l i reduce to 



A(x,z) ~ A(x) (l + 3e~ z/Hb ) (e zlHb - if 



(22) 




Fig. 1. Normalized density (solid) and Alfven velocity (dashed) 
as a function of the transverse coordinate x, which in dimen- 
sional units is given in Mm. 



AM 



vw 2 // 3 



-1/2 3/2 

96y AoPo Poo 



dpo 
dx 



(23) 



We see that, when z is a few times larger than H b , the wave 
amplitude is proportional to exp (-Ae 3z / Hb ), so that wave damp- 
in g occurs much faster than in the one-dimensional case studied 
by iHevvaerts & Priestl dl983b , with the characteristic damping 
length proport ional to In v. Qualitativ ely this result is in agree- 
ment with the lRuderman et all (Il998h solution, however the ex- 
act expressions for A and A that we obtained are different from 
their corresponding expression (see their Eqs. (45)-(46)). 

To obtain the equilibrium with constant Alfven speed we 
take H p -> ±tf b . Then Eqs. <[T9j»— ([20j reduce to 



A(x, z) ~ A(x) 



2z 



(e-^-3)(l-e-^) + ^ 



A(x) 



oj 2 vH? 



16V, 



APo 



dpo 
dx 



(24) 



(25) 



where Va(x) is no longer a function of z, as the divergence of 
the magnetic field exactly balances the effect of stratification. 
We see that, when z > H b , A oc z so that wave damping oc- 
curs much slower than in the IHevvaerts & Priestl (Q983) equi- 
libriu m where A oc z 3 . U sing the incorrect general solution for 
A. iRuderman et all ( 1998) obtained that, in the equilibrium with 
constant Alfven speed, the wave amplitude tends to a non-zero 
constant as z — > oo, whereas our corrected analytical solution 
shows that this is clearly not the case. See Appendix B for de- 
tailed derivations of these equilibrium solutions. 

In the following sections we will compare our numerical cal- 
culat ions to our corrected g eneral solution as well as to the orig- 
inal lRuderman et all dl998h general solution. This is done solely 
to emphasize the importance of the correction presented in this 
study. The authors note that, i n spite of an error in the calcula- 
tions of Ruderman et al .(1998), their results remain qualitatively 
correct. In particular, their conclusion that magnetic field diver- 
gence can strongly enhance phase mixing, is in agreement with 
the conclusions of this study. 



4. The Numerical Model 

We developed and tested a parallelized 2.5D visco-resistive lin- 
ear MHD code, consisting of a centered finite difference scheme 
(6th order in space), combined with a Runge-Kutta time step 
(4th order in time). The 2.5D nature of the code means that all 
variables are function of x and z only, however B and V have 
all three components. We then used this code to solve Eqs. (f3]l- 
©, written in dimensionless form. We used reference length, I, 
time, t, velocity, u, number density, N and mass density, R, to in- 
troduce the dimensionless quantities. For a typical coronal hole 
/ = 1 Mm, t = 1 s, u = 1 km s _1 , N = 1 irT 3 and R = 1 kg trT 3 . 
From here on, unless specified, all quantities are given in dimen- 
sionless form. 

In our numerical calculations we modeled various coronal 
structures, determined by Eqs. (TT3T> — (TTTTt, along with 



po(x) = p 



1 + tanh 



/x + 5 
"I 0.4 



- tanh 



(26) 



where po is the background density at x = 2, z — 0, which, using 
the mean ion mass (measured in proton mass units) p = 1.27 
and the electron number density «o = 6.43 x 10 14 , gives po w 
1.4 x 10~ 12 . The variations of po(x) and the Alfven speed across 
x at z = 0, are shown in Fig. Q] This figure shows a three fold 
increase in density, along with a corresponding 1 / V3 decrease 
in Alfven velocity from * 750 to » 350, as we move across the 
density boundary. This is in agr eement with observ ations of solar 
coronal plume boundaries (see Defore sTet alj|1997h . Fig. Q] also 
shows the density boundary has a characteristic scale of varia- 
tion, or half width, of xo ~ 1. For large magnetic scale heights, 
H\, > 700, Eq. ( f26b represents the boundary of a coronal plume. 
For lower magnetic scale heights, where H b < 100, we use it to 
represent the boundary of a divergent coronal loop. 

The density scale height was fixed at H p = 50 in all numer- 
ical calculations with stratified equilibria. In dimensional form 
H p = 50 Mm, which corres ponds to a temperature equal to 
2 MK (see Aschwanden 2004). The transverse density function, 
Eq. d26| i, was chosen such that the centers of the density bound- 
ary and magnetic field divergence coincide at x = 0. In Fig. 2 
the curvilinear coordinate system is shown, where the white box 
shows the size of the computation domain. 
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-10 -5 5 10 

Fig. 2. Magnetic surfaces defined by curvilinear coordinates 
if/(x,z) (dashed) and <p(x,z) (solid). The width of the computa- 
tional domain is shown by the white box. 

The plasma, which initially is at rest, is driven at z — 0, by 
B y = A Boo sin (cot) (l - exp [- (f/f ) 3 ]} , (27) 

V y = -A V A sin (cot) {l - exp [- (f/f ) 3 ]} , (28) 

where Boo is the magnitude of equilibrium magnetic field at 
x = z = 0, Va is the Alfven velocity at the driving boundary, 
Ao = 0. 1 is the dimension less driving ampl itude, t is the simu- 
lation time and fo = 100. iHood et alj d2002f) describe how the 
leading and trailing edges of a finite wave train develop into 
gaussian pulses that dissipate algebraically according to z~ 3 / 2 . 
Since these leading edge pulses dissipate at a much slower rate 
than the exponential rate we were seeking, it was important to 
minimize their effect on our numerical calculations. This was 
achieved by ramping up B y and V y over 20 wavefronts, by in- 
cluding the exponential dampening terms exp [— (f/fo) 3 ] seen in 
the boundary conditions of Eqs. (|27]>-(|28l. This minimizes the 
development of leading edge gaussian pulses and enabled the 
numerical calculations to reach a steady state before tracking of 
the wave amplitude begins. 

In all calculations the kinematic viscosity was v = 5 x 10" 5 
(anomalous), while the wave frequency was / = 0.1 (to — Q.2n). 
In dimensional form these correspond to v = 5 x 10 7 m 2 s _1 and 
/ = 0.1 Hz. See Sect. 6 for a discussion on the plausibility of 
these parameters. 

The size of the numerical domain was -2 < x < 2, 
< z < 100. We used a combination of zero-gradient 
(x = +2), line-tied (z - 0) and masked (z = 100) boundary con- 
ditions. Masking allows waves to propagate through the upper 
boundary (z = 100) without reflection. To use this method we 
extended the computational domain in the z-direction by 50% 
to create a masking region where all perturbations are exponen- 
tially damped. The corresponding numerical resolution, which 
includes the masking region, was 400 x 7500. The numerical 
convergence test was carried out by repeating numerical calcu- 
lations for the coronal structure with the most divergent mag- 
netic field (H\, = 5, H p = 50), with double numerical resolution; 
800 x 15 000. Each numerical calculation was run for a time 
t = 1000, on a dual-core Intel Xeon processor for approximately 
two weeks (depending on Hb), while the convergence check was 
run on 8 AMD Opteron processors for a period of one month. 
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Fig. 3. Variation of wavelength with height for (a) non-stratified, 
Hp = oo, and (b) stratified, H p = 50, weakly divergent coro- 
nal structures along x — 0. Solid, dotted, dashed and dashed- 
dotted lines correspond to Hb = 40, Hb = 100, Hb = 700, and 
Hb = oo respectively. In dimensional units, the wavelength, A\\, 
and height, z, are given in Mm. 



5. Numerical Results 

In this section we compare the analytical solutions of Sect. 3 to 
numerical calculations of Alfven waves propagating in divergent 
stratified coronal structures. We divided our numerical calcula- 
tions into two regimes: weakly divergent (Hb = 40, 100, 700, oo) 
and strongly (Hb = 5,10) divergent coronal structures. In both 
regimes, coronal structures both with and without density strati- 
fication (H p = 50) were considered. 

5. 1 . Variation of wavelength 

We begin by examining Fig. [3}a) and Fig. Ob), which show 
the variation of wavelength, /ty, with z for a range of magnetic 
scale heights, in both stratified and non-stratified weakly diver- 
gent coronal structures. The wavelength varies with height as 
the Alfven velocity, Eq. ( fTTb , depends on the local magnetic 
field strength and local plasma density, which were both chosen 
to decrease exponentially with height accor ding to their respec- 
tive sc ale heights, Hb and H p . In the standard lHevvaerts & Priestl 
(I1983I) phase mixing solution, Eq. d2Tl ). the Alfven velocity, and 
hence wavelength, is a function of x only; it does not change 
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Fig. 5. Alfven wave amplitude as a function of height along x — for non-stratified, H p = oo, coronal structures with magnetic field 
divergence (a) H\, = 40, ( b) H\, = 100, (c) H\, = 700, (d) H\, = oo along x = 0. The solid and dotted lines shows calculations from 
the corrected and original R uderman et alj d!998l) solution respectively, while the diamonds shows the numerical calculations. The 
dimensional units of amplitude, V y , and height, z, are respectively km s _I and Mm. 



with height. Whereas in our phase mixing solution given by 
Eqs. ||191l-ll20t, the Alfven velocity and wavelength are functions 
of both x and z. A divergent magnetic field will therefore cause 
the Alfven velocity, and thus the wavelength, to decrease with 
height. As a result stronger transverse gradients form, leading to 
greater viscous dissipation than would otherwise occur by stan- 
dard phase mixing alone. We therefore define enhanced phase 
mixing to mean phase mixing occurring in divergent magnetic 
fields which results in stronger dissipation, due to a reduction 
in the wavelength. In Fig. [4] we show an example simulation of 
the enhanced phase mixing effect for a stratified weakly diver- 
gent coronal structure, at times t = 200 and t = 600. Comparing 
Fig. HJa) and Fig. |4|b), we clearly see that as the Alfven wave 
propagates along the z-axis, the wavefront is strongly dissipated 
around x = 0, where the transverse density gradient is strongest. 

We note that the reduction in wave amplitude with height, 
due to enhanced phase mixing, results from two separate effects; 
firstly from the reduction in wavelength and the associated in- 
creased dissipation (as discussed above), and secondly from the 
area-divergence, or geometrical spreading, of the wave propagat- 
ing in a divergent magnetic field. We also note that the charac- 
teristic transverse scale does not approach zero, even in strongly 
divergent coronal structures; as the transverse scale gets smaller, 



stronger transverse gradients build up until a balance is eventu- 
ally reached with viscous dissipation. In the kinetic regime, the 
transverse scale can, as a result of strong phase mixing, be re- 
duced down to the ion-cyclotron radius. 

We see from Fig. Oa) that, in non-stratified coronal struc- 
tures, increasing the divergence of the magnetic field, by lower- 
ing H\,, increases the rate at which the wavelength shortens with 
height. When H\, = 40 the wavelength is reduced by approxi- 
mately a factor of 4 from A]\ ~ 5.5 at z = 0, to A\\ ~ 1.5 at z = 50. 
From Eq. ( fT~3b , where A oc to 2 oc 1 /A 1 , we see that a factor of 4 
reduction in wavelength increases the dissipation rate by a factor 
of 16. In the stratified coronal structures of Fig. [3fb) the picture 
is less straightforward; when H\, < 100 the wavelength again re- 
duced with height, but not to the same extent as in non-stratified 
coronal structures. When H\, = 40 the wavelength is reduced by 
approximately a factor of 2 from A|| w 5.5 at z — 0, to * 2.5 at 
z = 50. For H\, = 100 the wavelength remained unchanged, cor- 
responding to the special case where H p = \H\,\ see Eqs. d24"l) - 
( |25] |. This results from a constant Alfven velocity with height 
where the reduction in Alfven velocity, due to the divergence of 
the magnetic field, is exactly balanced by an increase in Alfven 
velocity due to density stratification. This specific configuration 
does not exhibit the wavelength reduction associated with en- 
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Fig. 6. Alfven wave amplitude as a function of height along x — for stratified, H p = 50, coronal structures of magnetic field 
divergence (a) = 40, ( b) H\, = 100, (c) H\, = 700, (d) H\, = oo along x = 0. The solid and dotted lines shows calculations from 
the corrected and original R uderman et alj d 1998b solution respectively, while the diamonds shows the numerical calculations. The 
dimensional units of amplitude, Vy, and height, z, are respectively km s _I and Mm. 



hanced phase mixing, but it does still exhibit the effects of area- 
divergence. When H\, > 100, the wavelength is seen to increase 
with height, due to the Alfven velocity's inverse dependence on 
density, Va 1 / \fpo, which exponentially reduces with height 
as a result of stratification. In Fig. Oa), H\, — oo is equivalent to 
a coronal structure perm eated by a uniform magne tic field, and 
thus represents standard He y vaerts & Priest! d 19831) phase mix- 
ing. Comparing this to Fig. [2b), where — oo represents phase 
mixing in a stratified coronal structure, we see that density strat- 
ification works to increase the wavelength, and thus reduce the 
overall dissipation rate. 



5.2. Weakly divergent coronal structures 

In Figs. [5] and [6] the numerically calculated Alfven wave ve- 
locity amplitude, V y , is plotted as a function of height for 
both non-stratified and stratified weakly divergent coronal struc- 
tures. These plots represent cross-sections of the wave's ampli- 
tude along x = 0, where the maximum dissipation occurs. We 
have also plotted our analytical solutions given by Eqs. (|19t - 
lEO), as well a s the c orresponding equilibrium solutions of 
Ruderma n et al.l ( 1998b. We h ave chosen to plot the incorrect 
original iRuderman et al.1 d 19981) solution to demonstrate the im- 



portance of our correction to the enhanced phase mixing solution 
of Sect. 3. Also note that we have only plotted V y and not B y , 
since the behaviors of these two quantities are identical for non- 
stratified structures, and differ only slightly in stratified struc- 
tures due to the amplification effect (see below). 

Firstly from Figs.[5]and[6] we note that lowering increases 
the rate of dissipation (see Sect. 5.1). Secondly we see that in 
stratified structures only, there is an amplification of V y at low 
heights, along with a corresponding reduction in B y (not shown). 
This well known effect results from the wave amplitude's de- 



pendance on density; Vy oc pJ /4 a nd B y oc p 1 ^ 1 * (see, e.g., 



,1/4 



IWright & Garm an 1998; Moran 2001). The reduction in density 
and subsequent increase in Alfven velocity associated with strat- 
ification, therefore leads to the amplification of V y seen in Fig. [6] 
We note that by using very high values of viscosity (v > 10~ 3 , 
corresponding to v > 10 9 m 2 s _1 ), it is possible to remove this 
wave amplification, as viscous dissipation would dominate any 
amplification at low heights. Thirdly we see that, for a fixed mag- 
netic scale height, Hb, Alfven waves dissipate slower in strati- 
fied structures than in non-stratified structures. This is due to the 
Alfven velocity decrea sing faster in non - strat ified structures than 
in stratified structures. IRuderman et al.l (1 19981) suggested that at 
low heights, harmonic Alfven waves dissipate according to the 
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Fig. 4. An example enhanced phase mixing simulation at times 
(a) t = 200 and (b) t = 600, for a weakly divergent, H b = 40, 
stratified, H p = 50, coronal structure. In dimensional units, x and 
Z are given in Mm, while V y is given in km s _1 . 



standard exp(-z 3 ) rate, since at these heights the reduction in 
wavelength required by enhanced phase mixing is not signifi- 
cant. This is confirmed by both Figs.[5]and[6]which show that, at 
low heights (z < 20), the waves dissipate independently of H\>. 
At larger heights (z > 20), the enhanced phase mixing mecha- 
nism begins to differentiate among the different equilibrium con- 
figurations. 

Next we compare our numerical calculations to our analyt- 
ical solutions given by Eqs. (fT9b — (l20l>, as well as to the cor- 
responding equilibrium solutions of Ru derman etafl (fl998h . In 
Figs. |3d) and |6fd), representing phase mixing in a uniform 
magnetic field, we see strong agreement between the numeri- 
cal model and both analytical solutions. As the magnetic scale 
height i s lowered from H h = o o to H b = 40, the analytical solu- 
tions of Ruderman et al. ( 1998) begin to diverge from the numer- 
ical model results, whereas our analytical solutions continue to 
show good agreement over a range of magnetic scale heights. 
This is most apparent in the stratified numerical calculations, 
Fig. [6] where the lower overall dissipation rates enable the differ- 
ence between standard and enhanced phase mixing to be more 
clearly seen. Notably when H\, = 40 for the stratified numerical 
model, Fig. |6|a), the Alfven wave is almost c ompletely dissi- 
pated by z ~ 60; half the height at which the iRuderman et al.l 
( 1998) analytical solution predicts and within one density scale 
height, Hp. 

5.3. Strongly divergent coronal structures 

We now consider coronal structures with strongly divergent 
magnetic fields, where the magnetic scale height is of the order 
of the characteristic transverse length scale, (xo/Hb ~ 1). This 
clearly violates the first assumption, xo/H «: 1, used to derive 
the general analytical solution (See Sect. 3). We would there- 
fore expect to see a large difference between the numerical and 
analytical solution. 

In Figs. [7] and [8] we plot the numerically calculated 
Alfven wave amplitude, V y , along x — for respectively, non- 
stratified and stratified strongly divergent coronal structures. As 
in Sect. 5.2 we also plotted our analytical solutions given by 
Eqs. (fT9l-(l20li. as well as the corresponding equilibrium ana- 




Fig. 7. Alfven wave amplitude as a function of height along 
x — for non-stratified, H p - oo, with strongly divergent 
magnetic field: (a) H\, = 5, (b) H\, = 10. The solid and 
dotte d lines shows calcu lations from the corrected and origi- 
nal l^dermanetal] (1 19981) solution respectively, while the dia- 
monds shows the numerical calculations. The dimensional units 
of amplitude, Vy, and height, z, are respectively km s and Mm. 



lytical solutions of Ruderma n et al.l (Il998l) . Comparing the two 
plots, we see that the effects of stratification are less pronounced 
for H\, = 5, than for H\> = 10. This occurs as the wave ampli- 
tude's dependence on magnetic scale height, exp(- exp{z/H\,)), 
dominates over the wave amplitude's dependance on viscosity 
and frequency, exp(-vw 2 ), in highly divergent magnetic fields. 
We therefore focus on Fig. [8] where once again we see that de- 
creasing Hb increases the dissipation rate as expected. We see 
that, when H\, = 5, we can fully dissipate Alfven waves with 
/ = 0.1 within z « 20, which is over six times low er than would 
occur in the standard iHevvaerts & Prie st (1983) phase mixing 
case; Fig. [5jd). This is also hal f the density scale heig ht H p , as 
well as half the height that the Ruder man et al.l d 19981) analyti- 
cal solutions predict, which as in the weakly divergent coronal 
structures, differs strongly from both our analytical and numeri- 
cal calculations. 

Note that the amplification of V y due to density stratification 
(see Sect. 5.2), is less pronounced than in the coronal structures 
with weakly divergent magnetic field seen in Fig. [6] This is be- 
cause the strongly divergent magnetic field very quickly gener- 
ates strong dissipation at heights much lower than in the weakly 
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Fig. 8. Alfven wave amplitude as a function of height along x = 
for stratified, H p = 50, with strongly divergent magnetic field: 
(a) H\, = 5, (b) //b = 10. The solid and d otted lines shows calcu- 
lations from the corrected and original Ikuderman et al.l (l998) 
solution respectively, while the diamonds shows the numerical 
calculations. The dimensional units of amplitude, V y , and height, 
Z, are respectively km s _1 and Mm. 



Fig. 9. Heating power along a flux tube for (a) non-stratified, 
Hp = oo, and (b) stratified, H p = 50, divergent coronal struc- 
tures. Solid, dotted, dashed, dot-dashed, triple-dot-dashed and 
long dashed lines correspond respectively to H b — 5, H\, = 10, 
H h = 40, H h = 100, H b = 700 and H b = oo. The dimensional 
units of heating power, E H , are xl0~ 18 J m~ 3 s _1 , while those for 
H b and z are Mm. 



divergent case, which is then able to overpower the stratification 
effect. Surprisingly the numerical results are still in good agree- 
ment with the corrected analytical solution along x = 0, even 
when the magnetic scale height is lowered to H b = 5 (xo/H b * 
0.2). It is interesting to note that standard phase mixing modeled 
by PIC (Particle in Cell) codes in the kineti c regime, gives dis- 
sipatio n rates in strong agreement with the |Hevvaerts & Priestl 
(119831) MHD dissipation rates (see iTsiklauri et al.l l2005h. The 
effectiv e resistivity seen in this paper, as well as in iTsiklauril 
(2006b), were also found to be many orders higher than the 
classical Braginskii value. This result along with the demonstra- 
tion that two physical descriptions, valid on completely different 
scales, can lead to the same dissipation rates, strongly supports 
the use of MHD and anomalous viscosity in this study. 

5.4. Viscous heating 

We now discuss the viscous heating generated by the enhanced 
phase mixing of Alfven waves in divergent coronal structures. 
This heating results from the phase mixing mechanism convert- 
ing the large scale wave energy to the small scale thermal en- 



ergy of the plasma, via viscous (or resistive) dissipation. In this 
section, we calculate the generated viscous heating power, En, 
and viscous heating flux, F^, from the numerical calculations of 
Sects. 5.2 and 5.3. 

The viscous heating power describes the location and mag- 
nitude of the energy dissipation, and is calculated from 



E H (x,z) = vp (VxV)V 



2 e zlH b 



(29) 



The viscous heating flux describes the flow of energy with the 
upward propagating waves, and is derived from the Poynting 
vector 



ExB 



(30) 



where subscripts and 1 refer respectively to equilibrium and 
perturbation quantities (except for yUo which is the magnetic per- 
meability), and E is the electric field given by E = -V\ x Bq. 
Substituting this into Eq. ( 1301 gives an alternative form for the 
wave Poynting vector 



B x x(Fi xfi ) 
Co 



(31) 
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Fig. 10. Heating flux along a flux tube for (a) non-stratified, 
H p = oo, and (b) stratified, H p = 50, divergent coronal struc- 
tures. Solid, dotted, dashed, dot-dashed, triple-dot-dashed and 
long dashed lines correspond respectively to Hb = 5, Hb — 10, 
H b = 40, H b = 100, H h = 700 and H b = oo. The dimensional 
units of heating flux, Fh, are xl0~ 12 J trT 2 s _1 , while those for 
Hb and z are Mm. 



from which we derive the wave Poynting flux, (5) = 
\B y Vy\Bo/fiQ, and finally heating flux, Fu = Fq - (S), as a func- 
tion of height 



\B v Vv\Bq ,„ 

F H (x,z) = Fo- y ' - e m \ 



(32) 



where Fq = AqB^Va/hq is trie wave Poynting flux at the coor- 
dinate origin, x = z — 0. Note that the heating rate and heating 
flux oscillate in time at twice the wave frequency, which ideally 
would be removed by time averaging Eqs. d29b and ( 1321 ). Instead, 
for numerical reasons, we have taken the wave envelope, scaled 
by 0.5, to be equivalent to time averaging when plotting Figs. [9] 
and[l0] The exp(z/ifb) terms seen in Eqs. d29l and d32l >. compen- 
sate for the area divergence of each flux tube, so that the effects 
of standard and enhanced phase mixing on wave dissipation, can 
be directly compared. 

In Figs. l9l and [TUl we show how the viscous heating power 
and viscous heating flux along a flux tube, depend on the mag- 
netic scale height. We begin by comparing Figs.|9|a) and |9jb), 
where we see that the effects of density stratification on the 
magnitude and location of peak viscous heating power, is mi- 



nor. There is a slight reduction in the magnitude of peak heat- 
ing power and lowering of the height at which this occurs, but 
these effects become increasingly insignificant as the magnetic 
field becomes more divergent {H\, is decreased). Therefore we 
focus on the more realistic case of Fig. 0b), where the effects 
of enhanced phase mixing are clearly apparent; increasing the 
magnetic field divergence, by reducing H b , causes; (i) the mag- 
nitude of heating power to increase, and (ii) the location of the 
heating power peak to be lowered. Indeed, comparing a strongly 
divergent magnetic field, Hb = 5, to a uniform magnetic field, 
H\, = oo, we see that the viscous heating power is increased by 
a factor of ten, from En « 4.2 x 10 13 to E H * 4.4 x 10 14 . This 
is approximately ten times the corresponding heating power re- 
quirement for an active region; Lh ~ 4.0 x 10 13 . We also see that 
the location of peak heating is lowered by a factor of three, from 
z*40toz^l3. 

Next we compare Figs. [TOl' a) and [TOl'b), where we see that 
density stratification affects weakly divergent coronal structures 
far more than their strongly divergent counterparts. This oc- 
curs because, in strongly divergent coronal structures, the wave- 
length reducing effects of enhanced phase mixing dominate over 
the wavelength increasing effects of density stratification (See 
Sect. 5.1). We again focus on the more realistic stratified case, 
Fig- QZHb), where firstly we note that the maximum viscous 
heating flux carried by the Alfven waves along a flux tube, is 
Fu =s 2.1 x 10 15 . Given that the heating flux requirement for 
an active region is Fu m 2.0 x 10 15 , this agrees well with our 
expectation that Alfven waves carry sufficient energy to heat 
the corona. Now comparing a strongly divergent magnetic field, 
Hb = 5, to a uniform magnetic field, Hb = oo, we see that the 
height at which the Alfven waves have dissipated 95% of their 
heating flux, which we take to be representative of the heating 
length scale, Lh, is lowered by a factor of six; from Lh ~ 90 
to Lh ~ 15. This follows from Sect. 5.3 where we showed that 
Alfven waves could be fully dissipated over six times lower in a 
strongly divergent field, than in a uniform magnetic field. 

5.5. Heating length scale 

In the previous section we demonstrated that for the parame- 
ters chosen in this study, which were based on observational evi- 
dence, Alfven waves contain sufficient energy to fulfill the coro- 
nal heating requirement. The only question is therefore at what 
height is this fulfillment achieved? To answer this, we consider 
the effects of altering the driving frequency, oj, and the shear 
viscosity, v, on the heating length scale, Lh, generated by the en- 
hanced phase mixing of Alfven waves propagating in divergent 
and stratified, H p = 50, coronal structures of varying magnetic 
field divergence. We use Eqs. ( fT9l and ( l20l i. from the analytical 
solutions of Sect. 3, which were found to be in strong agreement 
with the numerical calculations of Sects. 5.2 and 5.3, to perform 
a parametric study of the variation in heating length scale with 
driving frequency and shear viscosity. In the previous section the 
heating length scale, Lh, was defined to be the height at which 
95% of the Alfven wave Poynting flux, along a flux tube, had 
been dissipated. For an active region, this heating length scale 
is also required to be within a density scale height, which in 
ou r model corresponds to L h < 50, to agree with an assertion 
bv lAschwanden etakl (|2000) that the heating scale height is less 
than the density scale height. 

In Fig. [TT] we plot the heating length scale as a function of 
o> and v, for magnetic scale heights, Hb = 5, 10, 40, oo. Firstly, 
from the standard phase mixing case shown in Fig. fTTT d). we see 
that to generate the required heating length scale of Lh < 50 us- 
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Fig. 11. Heating length scale as a function of driving frequency, co, and shear viscosity, v, where (a) H\, - 5, (b) H\, - 10, (c) 
H\, = 40 and (d) = oo. Solid, dotted, dashed and dot-dashed lines correspond to Lh = 25, Lh = 50, Lh = 75 and Lh = 100 
respectively (H p = 50). The dimensional units of Lh and H\, are Mm, co are radians s~ l and v are xlO 12 m 2 s . 



ing observable, co w 0.01, Alfven waves, requires an anomalous 
viscosity 8 orders of magnitude higher than classical Braginskii 
viscosity, v w 2 x 10~ 9 . This demonstrates the main problem with 
standard phase mixing in the corona; we need to use extremely 
high values of shear viscosity, to dissipate Alfven waves low 
enough in the corona to contribute to coronal heating. Secondly, 
we see that even for high frequency, co > 1 .0, Alfven waves, we 
still require the viscosity to be 3 to 5 orders of magnitude higher 
than the Braginskii value. Clearly, without the combination of 
anomalous viscosity and high frequency Alfven waves, standard 
phase mixing is not a viable coronal heating mechanism. 

In Fig. [TTT c). we see that the enhanced phase mixing of 
Alfven waves in weakly divergent magnetic fields, = 40, 
does not significantly alter the heating length scale. Whereas 
in Figs. [TTT a) and [TTT b). which represent the enhanced phase 
mixing of Alfven waves in strongly divergent magnetic fields, 
H\, < 10, the heating length scale is significantly reduced. For 
H\, = 5, we see that the enhanced phase mixing of observ- 
able, co w 0.01, Alfven waves generates the required heating 
length scales, Lh ^ 50, using classical Braginskii viscosity. 
As in Sect. 5.3, this occurs because in highly divergent mag- 
netic fields, the wave amplitude dependence on magnetic scale 
height, exp(- exp(z///b)), dominates over the wave amplitude 
dependance on viscosity and frequency, exp(-va) 2 ). Comparing 



Figs. [TTT a) and [TTTd). we see that in standard phase mixing, 
small changes of co or v have large effects on the heating length 
scale, whereas the opposite is true for enhanced phase mixing in 
a strongly divergent magnetic field. We also see that, for a given 
wave frequency, to generate a heating length scale of Lh < 50 
using enhanced phase mixing, requires a shear viscosity eight 
orders of magnitude lower than standard phase mixing. Finally, 
for co ~ 0.6 and v = 5 x 10~ 5 , corresponding to the values used 
in this study, the magnetic scale height must be H\, < 40 for the 
heating length scale to be Lh < 50, which is in agreement with 
the corresponding numerical calculations of Sect. 5.4. 

Note that we must exercise caution when considering low 
frequency Alfven waves, since the derivation of the analytical 
solution of Sect. 3 assumes the wavelength is smaller than or of 
order the transverse scale, xq. To demonstrate the accuracy of 
Fig. Q~T] we therefore conducted numerical calculations for the 
low frequency case, where co ~ 0.01 and v = 2x 10~ 9 for H b = 5, 
H p = 50. We found the calculated heating length scale, Lh ~ 50, 
to be in good agreement with the results seen in Fig.fTTl 



6. Discussion and Conclusions 

We analytically and numerically studied the phase mixing 
of Alfven waves propagating in weakly divergent, H\, = 
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40, 100,700,oo Mm, and strongly divergent, H b - 5, 10 Mm, 
stratified, H p = 50 Mm, coronal structures. These numerical cal- 
culations were used to validate our analytical solution, which 
was obta ined by correcting an er ror in the general analytical so- 
lution of iRuderman et all dl998l) . For convenience, in this sec- 
tion we revert to using dimensional units. 

We began in Sect. 5.1 by showing that density stratification 
and magnetic field divergence are two opposing factors affecting 
the wavelength; stratification works to increase the Alfven veloc- 
ity while divergence works to decrease it. An increasing Alfven 
velocity increases the wavelength leading to larger transverse 
scales, and thus reduces the wave dissipation rates. Conversely 
a decreasing Alfven velocity shortens the wavelength leading to 
smaller transverse scales, and thus increases the wave dissipa- 
tion rates. Note that in both cases, we assume that the angular 
wave frequency, a>, remains fixed. We used this to define the 
concept of enhanced phase mixing as; phase mixing occurring 
in divergent magnetic fields, which results in stronger dissipation 
due to a reduction in the wavelength. We found that in stratified 
coronal structures, enhanced phase mixing occurs only when the 
magnetic scale height is less than twice the density scale height, 
Hb/H p < 2. Therefore the enhanced phase mixing of Alfven 
waves in the corona will only occur when Hb < 100 Mm, given 
that the density scale height is H p = 50 Mm. This also means 
that enhanced phase mixing does not occur in typical coronal 
plumes where H p «: Hb, but is instead limited to the density 
boundaries of highly divergent coronal structures, e.g. coronal 
loops and arcades. 

In Sect. 5.2 we compared our numerical calculations to 
our corrected analytical solution, Eqs. ([T0li-([T3j>. a s well as to 
the previous analytical solution of Ruderman et al. ( 1998). For 
very weakly divergent magnetic fields, corresponding to coro- 
nal plumes where Hb > 700 Mm, our analytica l and numerical 
calcul ations approximately agree with those of Ruderman et alj 
(1998f). Their conclusion that enhanced phase mixing is not 
an important dissipation mechanism in coronal plumes (a very 
weakly divergent magnetic field), is valid. However due to a 
derivation error ( "see S ect. 3), the dissipation rates given by the 
Ruderma n et alj d!998l) analytical solution show very small vari- 
ation with Hb, implying that magnetic field divergence has little 
effect on the dissipation of Alfven waves. The numerical calcu- 
lations we conducted in conjunction with our analytical solution 
show conclusively that this is not the case. Indeed for weakly 
divergent, Hb = 40 Mm, stratified, H p = 50 Mm, coronal struc- 
tures, we saw that 0. 1 Hz Alfven waves could be fully dissipated 
within one density scale height which is half the height that the 
original lRuderman et al.l d 19981) analytical solution predicts. 

In Sect. 5.3 we found that, the corrected analytical solu- 
tion, Eqs. (fTOb — (fT3b. can accurately describe highly divergent 
coronal structures, even when the thin structure approximation, 
that the ratio of the characteristic scales in x and z-directions is 
small (xq/H «: 1), used in its derivation begins to be violated. 
Importantly we found that 0.1 Hz Alfven waves propagating in 
strongly divergent, Hb — 5 Mm, stratified coronal structures can 
be full y dissipated within 20 M m, which is again half the height 
that the Ruderman et al. ( 1998) analytical solutions predict. This 
is als o over six times lower th an would occur as a result of stan- 
dard iHevvaerts & Pries! d 19831) phase mixing in uniform mag- 
netic fields and is less than half of the density scale height. This 
result suggests that the importance of enhanced phase mixing as 
a mechanism for dissipating Alfven waves in the solar corona (a 
stratified and divergent medium), has been seriously underesti- 
mated. 



In Sect. 5.4 we investigated the viscous heating that results 
from enhanced phase mixing. We found that 0.1 Hz Alfven 
waves of amplitude V y ~ 55 km s _1 propagating in strongly di- 
vergent, Hb — 5 Mm, stratified coronal structures, can generate 
viscous heating fluxes Fu w 2100 J irT 2 s , with the associated 
viscous heating power peaking at En w 4.4 x 1 0~ 4 J m~ 3 s~ 1 . This 
compares well to an active region's Fu ~ 100 - 2000 J irT 2 s _1 
and Eh ~ 4.0 x 10~ 5 J irT 3 s _I heating requirement (see 
Asch wandenl 12004). Crucially, we also found that, as well as 
increasing the magnitude of the heating power by a factor of 
ten, enhanced phase mixing lowers the heating length scale by a 
factor of six, when compared to standard phase mixing in a uni- 
form magnetic field. This is perhaps the most important effect, 
as it enables the Alfven waves to dissipate their energy at low 
heights in the corona, within a density scale height. We there- 
fore conclude that, given strongly divergent coronal structures, 
the enhanced phase mixing mechanism is a viable method for 
heating the corona to temperatures in excess of 1 MK. 

Finally in Sect. 5.5, we considered the dependence of the 
heating length scale, Lh, on the driving frequency, a), and shear 
viscosity, v. Overall we found that the enhanced phase mix- 
ing of Alfven waves in strongly divergent magnetic fields, can 
dissipate the wave energy within a density scale height, us- 
ing a shear viscosity up to eight orders of magnitude lower 
than required by standard phase mixing in a uniform magnetic 
field. Specifically, we found that for strongly divergent magnetic 
fields, Hb — 5 Mm, the enhanced phase mixing of observable 
u> = 0.01 rads s _1 Alfven waves, results in a heating length 
scale of Lh < 50 Mm, using classical Braginskii viscosity; we 
therefore do not need to invoke anomalous viscosity to dissipate 
observable Alfven waves in the corona. Therefore we conclude 
that the value of shear viscosity required for the enhanced phase 
mixing mechanism to be a viable method for coronal heating, 
strongly depends on the heating length scale of an active region. 

Our decision to focus on active regions was led by the claim 
bv lAschwanden etal] (120071) that 83% of the coronal heating re- 
quirement is located there. While our focus on strongly divergent 
coronal structures was necessary to meet the requirement that 
the heating length scale be less than the active region density 
scale height, Lh < H p . Given this, the potential of the enhanced 
phase mixing mechanism to generate significant heating to ful- 
fill the coronal heating requirement is dependent on three critical 
parameters: the wave amplitude, wave frequency and the value 
of shear viscosity. Our use of wave amplitude V v » 55 km s' 1 is 
strongly supported by o bservations (lMoranl200ltlBaneriee et ail 
1998: lDovle etal.ll9 98). Our decision to use an anomalous kine- 
matic viscosity value of v = 5 x 10 7 m 2 s _1 , whi ch is four 
orders of magnitude higher than that given by the Braginskii 
(1 19651 ) shear viscosity t ensor, is supported by numerical evi- 
dence, Tsiklauri (2006b), which suggests that this is indeed the 
effective value in the corona. We also chose to use a relatively 
high Alfven frequency, / = 0.1 Hz. Given the turbulent na- 
ture of the photosphere it would seem likely that high frequency 
Alfven waves, / > 0.1 Hz, do indeed exist; it has been sug- 
gested the y are generated by micro-fl aring at the chromospheric 
level (see lAxford & McKenzid 19921) . As to whether their power 
spectrum is sufficient to account for coronal heating is an open 
question. Estimations of the high frequency Alfven wave spec- 
trum can be made by projecting back f rom the observed spec- 
trum at 0.3 AU (see ICranmer et al]| 19991) . Observations of high 
frequency Alfven waves are extremely difficult, as the lack of in- 
situ measurements means it is not possible to measure directly 
the magnetic field component. Indeed in Sects. 5.2 and 5.3, we 
plotted the variation of the V y with height, since it is the com- 
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panion B y component that is actually measured, via doppler line 
broadening. Therefore observational studies of the variation of 
the Alfven wave B y component with height should, with suffi- 
cient temporal and spatial resolution, be able to confirm or re- 
fute our model of coronal heating via the enhanced phase mix- 
ing mechanism. The upcoming ROSA (Rapid Oscillations in the 
Solar Atmosphere) ground based instrument, will make simulta- 
neous observations of the photosphere, chromosphere, transition 
region and corona at high cadence, and may therefore be able to 
directly observe the propagation and dissipation of these Alfven 
waves. 

In this study, we have shown that the enhanced phase mixing 
mechanism can fulfill the corona heating requirement in a suf- 
ficiently divergent active region. In reaching this conclusion we 
note that our numerical model is relatively simple. Therefore the 
inclusion of additional physical effects which are known to en- 
hance the wave dissipation, such as pressure, non-linearality and 
3 -dimensionality, would result in further heat deposition from 
the driven Alfven waves. The heating results shown in this pa- 
per should therefore be considered as a lower limit. It is there- 
fore essential to conduct realistic numerical simulations of wave 
dissipation using observational, 3D magnetic fields with waves 
driven by observed doppler shift data. Indeed, simulations incor- 
po rating observe d magnetic fields have already been attempted 
by lOfmanl d2007l) . in relation to coronal seismology. The inclu- 
sion of physical realism into numerical simulations is crucial to 
fully understanding coronal heating, and must therefore form the 
basis of future work. 
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Appendix A: General solution 

Here we derive our corrected general solution, Eqs. (fTOb— dT~3l> , 
for the enhanced phase mixing of Alfven waves in weakly di- 
vergent stratifie d coron al structures. As the derivation of the 
Ruderma n et al.l 0998) general solution is valid up to their 
Eq. (36), we will not present the full derivation here. For clar- 
ity we begin at Eq. (33) of their analytical paper, 

dx 

Tup 



1 dip 
'l~d~z 



dz 



1 dip 
J dx 



(A.l) 

Now given that Bo = \Bq\ and assuming Bo ~ #0z = Boodip/dx 
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Bn 



dz 



and therefore 

~ B QO dz' H * Boo \ d<f> 
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iRudermanetafld 19981) 
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ti) yfcr 
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deb, 



We now make another approximation that d<p/dz » d<p/dx. The 
total derivative therefore becomes 



(d<f>\ (d<f>\ dcf> 
d(p = [dx- ) dX + Vd-z) dZKS Tz d ^ 



which can then be substituted into Eq. ( IA.5 

~ ue \ Ti ) ,\ dz '- 
J V A (x,z') 
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Now substitute this a long with Eq. iA.3i into Eq. (30) of 
iRudermanetakld 19981) 



to give 
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Now since 
d dxdV A d 



dtp dif/ dx 8V A 
Eq. dA.lOb becomes 
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Next we substitute this into Eq. ( |A.91 > to give 

i \ 2 



A (<D, iff) 



CO 



2y B 



I* 



dx dV A , 
-dz 



'I dip dx 



dz 
d(p' 



(A.13) 



Now since A = J A i//jd<b and using Eq. iA.6i once more, we 



arrive at 



A(x,z) = 



2B 00 



l J v A J v\ 



dx dV A 
f\ dip dx 



dz 



dz'. 



(A. 14) 



V0 

which is equivalent to our general solution, Eqs. (TTOb— (TT~3b. 



(A. 4) Appendix B: Equilibrium solutions 

The error in the derivation of the iRuderman et al.l (Q998) gen- 
eral solution also carries through into their equilibrium analyt- 
ical solutions. We therefore present here derivations for these 
corrected equilibrium solutions. In each case we begin from our 
(A. 5) general solution for the enhanced phase mixing of Alfven waves 
in weakly divergent stratified coronal structures, Eqs. (TTOt — (fT~3l>. 
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B.1. Non-divergent, non-stratified (H b - H p - oo) 



Bo = Boo, Po (x) = po O) . 
Substitute these equations into Eq. ( fTTT i to give 
Boo 



(B.l) 



Va'oPo 



Now substitute this into Eqs. (fTOb — <fT3T>; our general solution 

Mx,z) 
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where the integral / is given by 
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2BooVPo 0x 
Substitute this into Eq. ( 1B.3I ) 
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and then simplify to give 
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Since Bo = Boo, this can be given in terms of Va using Eq. ( IB. 2b 
co 2 v ldV A " 2 
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6Vl \ dx 



z 3 . 
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which is equivalent to Eq. (12Tb . 



B.2. Divergent, non-stratified (H b + H p = oo) 



Substitute this into Eq. ( IB. 10b 



Here we derive Eq. (|2 Tb: the solution first derived by 
iHevvaerts & PriesH (1 1983b for the standard phase mixing of 
Alfven waves in non-stratified coronal structures permeated by J ^ JC ' & 
uniform magnetic fields. From Eqs. (fl3Tl-(fT7Ti we see that 
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x Hb)\ dx 

which is equivalent to Eqs. (l22b— (f23l) . 
B.3. Divergent, stratified, (H b = \H p ) 

Next we derive Eqs. (f24b — d25b ; the solution for the enhanced 
phase mixing of Alfven waves in stratified coronal structures 
permeated by weakly divergent coronal magnetic fields where 
H b = \H p . From Eqs. (fl"5ll-([T7l> we see that 



(B-5) Bo(x,z) = Booe- zJH \ po (x, z) = po M e~ 2z/Hb . 
Substitute these equations into Eq. dTTb to give 
Boo 
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Note that this is a function of x only. Now substitute this into our 
general solution given by Eqs. (fTOb— dT3b 
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where the integral / is given by 

H h y/jJS I x \ dpo 
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Now we derive Eqs. (F22l -d23b: the solution for the enhanced 
phase mixing of Alfven waves in non-stratified coronal struc- 
tures permeated by weakly divergent magnetic fields. From Substitute this into Eq. dB~T6b 
Eqs. dT3J — dTvJ we see that 
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which is equivalent to Eqs. (f24b — (t25b. 
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B.4. Divergent, stratified (H b + \H P ) 

Finally we derive Eqs. (fl~9T> — <f20b; the solution for the enhanced 
phase mixing of Alfven waves in stratified coronal structures 
permeated by weakly divergent magnetic fields where H b + 
\H p . From Eqs. (IT5T>— (ITTb we see that 

Bo (x, z) = B {yQ e- zlH \ po (x, z) = po (*) (B.20) 
Substitute these equations into Eq. (flTT i to give 



V A (x,z) 



fi 0» r z(V2H„-l/H k ) 



-yWo 

Now substitute this into Eqs. (fTOt— (IT3T> 
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where the integral / is given by 
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Substitute this into Eq. ( IB. 22b 
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